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We present time-series photometry of stars located in the extremely young open 
cluster Berkeley 59. Using the 1.04 m telescope at ARIES, Nainital, we have identified 
42 variables in a field of ~ 13 x 13 around the cluster. The probable members of the 
cluster are identified using (V, V — I) colour-magnitude diagram and (J — H, H — K) 
colour-colour diagram. Thirty one variables are found to be pre-main sequence stars 
associated with the cluster. The ages and masses of pre-main sequence stars are derived 
from colour-magnitude diagram by fitting theoretical models to the observed data 
points. The ages of the majority of the probable pre-main sequence variable candidates 
range from 1 to 5 Myrs. The masses of these pre-main sequence variable stars are found 
to be in the range of ~0.3 to ~3.5 Mq and these could be T Tauri stars. The present 
statistics reveal that about 90% T Tauri stars have periods < 15 days. The classical 
T Tauri stars are found to have larger amplitude in comparison to the weak line T 
Tauri stars. There is an indication that the amplitude decreases with increase of the 
mass, which could be due to the dispersal of disk of relatively massive stars. 
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1 INTRODUCTION 

Variable stars are useful tools to improve our understand- 
ing on stellar evolution and structure of stars. Star clus- 
ters are unique laboratories to study the stellar evolution 
as star clusters provide a sample of stars having same age, 
distance, initial composition and spanning range of masses. 
Additionally, the light curves of several stars can be pro- 
duced simultaneously to know the real variable stars among 
them. Studies of variable stars in star clusters are carried 
out since several years (Herbst et al. 1994; Edwards et al. 
1993; Oliveira & Casali 2008; Grankin et al. 2008 ). Star 
clusters have different types of variable stars as stars show 
variability at various stages of their evolutionary phases. 

In the present study, we made a search for variable stars 
in an extremely young open cluster Berkeley 59 (Be 59). Be 
59 is a young cluster and it contains a significant numbers of 
pre-main sequence (PMS) stars i.e, stars evolving from their 
birthline to the main-sequence (MS). This cluster (02000 = 
OO' 1 02 m 13 s , (5 2 ooo = +67° 25' ll" , 1 = 118°.22, b = 5°.00) is 
located at the centre of OB4 stellar association surrounded 
by HII region, Sharpless region S171 (Yang & Fukui 1992). 
Previous studies have revealed that Be 59 consists of many 
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PMS stars having age around 2 Myr (Pandey et al. 2008). Be 
59 was first studied photographically by Blanco & Williams 
(1959). They concluded that stars of spectral types 07 to 
B5 are still surrounded by a part of the parent molecular 
cloud. A multiwavelength study on Be 59 was presented 
by Pandey et al. (2008, hereafter P08). They found that 
this cluster is located at a distance of 1.0 kpc and having 
a variable extinction of E(B — V) = 1.4-1.8 mag. They also 
found that around 32% of the Ha emission stars in Be 59 ex- 
hibit NIR excess indicating that they still have their inner 
disk. During a campaign of photometric monitoring, Ma- 
jaess et al. (2008) presented the nature of the three eclips- 
ing systems (BD+66°1673, 2MASS 00104558+6127556 and 
2MSASS 19064659+4401458) in the field of Be 59 using BV 
photometry and spectroscopic observations. From the above 
studies it is clear that Be 59 contains high, intermediate and 
low mass PMS stars. 

PMS objects are, in general, classified into T Tauri stars 
(TTSs) and Herbig Ae/Be stars. The TTSs are newly formed 
low- mass stars (<3M©) which are contracting towards the 
MS (Herbig 1977), while stars of intermediate mass (>3-10 
Mq) are known as Herbig Ae/Be stars (Herbig 1960; Finken- 
zeller & Mundt 1984; Strom et al. 1972). Herbig Ae/Be stars 
are either contracting towards MS or they have reached 
the MS. The TTSs can be classified further as Weak line 
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TTSs (WTTSs) and Classical TTSs (CTTSs) on the ba- 
sis of the strength of the Ha emission lines (Strom et al. 
1989). The strength of the emission line is measured by its 
equivalent width (EW). The WTTSs exhibit a weak, narrow 
Ha (EWSJ 10A) in emission with no or little infrared ex- 
cess, while CTTSs generally display strong Ha emission line 
(EW>10A), large ultraviolet and infrared excess. The large 
Ha emission (EW>10A) is believed to be a result of the in- 
creased mass accretion onto the star (e.g. Cabrit et al. 1990; 
Calvet & Hartmann 1992). It is found that PMS stars show 
strong variability both in photometric and spectroscopic 
(van den Ancker et al. 1998; Catala et al. 1999) observations. 
This could be due to the presence of circumstellar obscura- 
tion and spots present on the photosphere. In fact variabil- 
ity is a defining characteristics of TTSs (Herbst et al. 1994; 
Scholz et al. 2009 and references therein). Both WTTSs and 
CTTSs show variation in their brightness. These variations 
are found to occur at all wavelengths, from X-ray to in- 
frared. Variability time scale of TTSs ranges from few min- 
utes to years (Appenzeller & Mundt 1989). The photometric 
variations are believed to originate from several mechanisms 
like circumstellar gas and dust (remnant of parent molecu- 
lar cloud), accretion and magnetic field (Herbst et al. 1994). 
The variations in the brightness of TTSs are most probably 
due to the presence of spots (cool or hot) at stellar surface 
and circumstellar disk. 

The cool spots on the surface of the stars are produced 
by the emergence of stellar magnetic fields on the photo- 
sphere, and is thus an indicator of magnetic activity. As 
these spots are present on the photosphere, they rotate with 
the star. If these spots are symmetrically distributed over the 
photosphere, periodic variations are seen in the light curves 
of the stars. Spot configurations remain constant over time- 
scales of several days, but undergo changes on time-scales 
of weeks (Hussain 2002) . The cool spots on the photosphere 
of star are responsible for brightness variation in WTTSs, 
and these objects are found to be fast rotators as they have 
either thin or no circumstellar disk. Therefore the variability 
in these objects may be mainly due to the brightness varia- 
tions on the photosphere modulated by the stellar rotation 
(Messina et al. 2010). 

The hot spots on the surface of young stars are supposed 
to be formed on the surface by infalling gas and are thus a di- 
rect consequence of accretion (Lynden-Bell & Pringle 1974; 
Koenigl 1991; Shu et al. 1994). Irregular or non periodic 
variations are produced because of changes in the accretion 
rate. The time scales of varying brightness range from hours 
to years. The hot spots cover a smaller fraction of the stel- 
lar surface but higher temperature causes larger brightness 
variations (Carpenter et al. 2001). Thus, accreting CTTSs 
which are surrounded by circumstellar accretion disk and 
possess hot spots produced by accreting material from disk 
on to the star show complex behaviour in their optical and 
near infrared (NIR) light curves (Scholz et al. 2009). 

The Herbig Ae/Be stars also show variability as they 
move across the instability region in the Hertzsprung-Russell 
(HR) diagram on their way to the MS. The right combina- 
tion of mass, temperature and luminosity of these stars make 
them to pulsate. Several PMS pulsators have been detected 
until now (Breger 1972; Kurtz & Marang 1995; Marconi et 
al. 2000; Donati et al. 1997; Zwintz et al. 2005, 2009; Zwintz 
& Weiss 2006). Catala (2003) presented detailed review on 



Table 1. Log of optical photometric CCD observations. 'Exp.' 
and 'n' refer to exposure time and number of frames respectively. 
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the PMS pulsations and expected that PMS stars could pul- 
sate as S Scuti stars, however it is necessary to confirm their 
PMS nature because these stars are found above the MS 
where it is very difficult to say whether they are PMS or 
post MS stars. But in cluster environment where the age of 
the cluster is constrained well, the stars found towards the 
cluster region could be PMS stars. 

A search of photometric variables in young open cluster 
Be 59, which is found to have a rich population of PMS stars 
with a large mass range could give us a better insight of the 
variability time-scales, the amplitude of brightness variation 
etc. With this aim, we observed the cluster Be 59 at different 
epochs from 2006 to 2010. The paper is organised in the 
following manner. In Section 2, we describe the observations 
and data reduction techniques. We discuss membership of 
variable stars on the basis of photometric observations in 
Section 3. In Section 4, we determine periods of variable 
stars. We discuss variability characteristics as well as period 
and amplitude distribution of variable stars in Section 5, 
while in Section 6 we conclude our studies. 



2 OBSERVATIONS AND DATA REDUCTION 

Our observing campaign of Be 59 region covers 18 nights 
from 2006, December 21 to 2010, October 31. Fig. 1 shows 
the observed field of Be 59. The observations in V band 
were obtained using 2kx2k CCD camera attached to the 
104-cm ARIES telescope. As each pixel covers 0.38 arcsec 
over sky, the field of view is ~ 13 arcmin on each side. On 
each night, at least two frames of target field in V band were 
secured to cover long period variables. The observations of 
Be 59 consist of total 257 CCD images. In addition, bias 
and twilight frames were also taken on each night. In order 
to improve the signal to noise ratio the observations were 
taken in 2x2 pixel binning mode. Log of the observations is 
given in Table 1. The preprocessing of the CCD images was 
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Figure 1. The observed CCD image of Be 59 cluster region. The squares show the location of candidate variables identified in the 
present work. 



done by IRAF which includes bias subtraction, flat fielding 
and cosmic ray removal. The instrumental magnitude of 
each star was obtained by DAOPHOT software package 
given by Stetson (1987). To get the magnitude of each 
star both the aperture and PSF photometry have been 
done because the PSF photometry yields better results in 
crowded region. The standardisation of the observations 
was done by observing the standard field SA 98 (Landolt 
1992). Instrumental magnitudes were converted to standard 
magnitudes using the following transformations equations 



V + (4.448 ± 0.004) - (0.037 ± 0.006)(V - /) + 0.15 * Q 



i = I+ (4.767 ± 0.008) - (0.04 ± 0.005)(V - I) + 0.05 * Q 

where v and i are the instrumental magnitudes and Q is the 
airmass. 



2.1 Comparison with Previous Photometry 

A comparison of the present photometry with that of P08 
yields 342 common stars. Fig. 2 plots difference A (present 
data— literature data) as a function of V magnitude. The 
comparison indicates a decreasing trend in AV values. The 
AV is w —0.13 mag upto V«14 mag which decreases to w 



—0.08 mag upto V^IS mag. Whereas present (V— /) colours 
are comparable to those given by P08. 



2.2 Identification of Variables 

The DAOMATCH (Stetson 1992) routine of DAOPHOT 
was used to find the translation, rotation and scaling solu- 
tions between different photometry files, whereas DAOMAS- 
TER (Stetson 1992) matches the point sources. In order to 
remove frame-to-frame flux variation of stars due to airmass 
and exposure time, we used DAOMASTER programme to 
get the corrected magnitude. This task makes the mean flux 
level of each frame equal to the reference frame by an ad- 
ditive constant. Present observations have 257 frames and 
each frame corresponds to one photometry file. DAOMAS- 
TER cross identified 472 stars in different photometry files 
and listed their corrected magnitudes in .cor file, which was 
used for subsequent tasks. The mean value of the magnitude 
and square root of the variance (RMS) of the data were es- 
timated using the observations of each star. The RMS as a 
function of instrumental V magnitude (Vi ns t) is shown in 
Fig. 3, which indicates that the majority of the stars follow 
an expected trend i.e., S/N ratio decreases as stars become 
fainter. However, a few stars do not follow the normal trend 
and exhibit relatively large scatter. These could be either 
due to the large photometric errors or due to the variable 
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Figure 2. Comparison of the present photometry and photome- 
try given by Pandey et al. (2008). 
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Figure 4. The RMS along with mean photometric errors of the 
detected variable stars as a function of instrumental magnitude. 
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sample of 42 variables. Their identification numbers, coordi- 
nates and photometric data are given in Table 2. The RMS 
along with mean photometric errors of detected variables as 
a function of instrumental magnitude are shown in Fig. 4. 
The photometric error is found to be ~ 0.010 mag at Vi„ s t 
< 14 mag, whereas the value of photometric error increases 
to ~0.08 mag at ~18 mag. 

Out of three variable stars reported by Majaess et al. 
(2008), only one variable star (BD+66°1673) is located in 
the present observed field of Be 59. Unfortunately, this star 
was saturated in our observations. 

In order to find out the possible causes of photomet- 
ric variations in the identified variable stars in the region 
of Be 59, it is important to ascertain their membership. 
In the absence of spectroscopic observations, we established 
photometric membership on the basis of (V, V — I) colour- 
magnitude and (J — H, H — K) colour-colour diagram in next 
Section. 



Figure 3. RMS for each star as a function of brightness. Open cir- 
cles represent candidate variables identified in the present work. 
Some stars showing large RMS are not considered variable stars 
as they are located near the edge of the CCD. 



nature of the stars. We considered a star as variable if its 
RMS is greater than 3 times of the mean RMS of that bin. 
Forty five candidate variables were thus identified on the 
basis of above mentioned criterion. On the basis of a careful 
inspection of light curves and the location of the sources on 
the CCD frames, we selected a sample of 37 stars to study 
their nature. A few stars were rejected as they were lying 
near the edge at the CCD. The mentioned technique, how- 
ever, could not detect low amplitude variables and 5 vari- 
ables namely V14, V18, V21, V23 and V19 were detected 
manually by inspecting their light curves. Thus we have a 



3 CLUSTER MEMBERSHIP 

The colour-magnitude diagram (CMD) of the cluster region 
clearly reveals contamination due to field star population (cf. 
figure 7 by P08). P08 found that the cluster is located at a 
distance of 1.00±0.05 kpc and has a differential reddening 
E(B — V)=1.4-1.8 mag. They also identified a foreground 
field star population at a distance of ~470 pc with a E(B — 
V)=QA0 mag (cf. figure 4 of P08). Therefore, the sample of 
of identified variables may be contaminated by foreground 
population. We have used V, (V - J) CMD and NIR J - 
H,H — K colour-colour diagram to find out the association 
of the identified variables with the cluster Be 59. 
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Table 2. Basic parameters of 42 variables in open cluster Be 59. Ho sources are taken from Pandey et al. (2008). 
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3.1 V , (V — I) Colour-Magnitude Diagram 

Fig. 5 displays V, (V — I) CMD for the stars observed in 
the cluster. The identified variables are marked by triangles 
(CTTSs) and filled circles (WTTSs) (cf. section 3.2). The 
zero-age main-sequence (ZAMS) by Girardi et al. (2002) cor- 
rected for the distance of 1.0 kpc and E(B — V) =1.45 mag 
is shown by thick continuous line. The pre-main-sequence 
(PMS) isochrones for 0.1, 1, 5 and 10 Myrs by Siess et al. 
(2000) have also been overplotted. The ZAMS corrected for 
a distance of 470 pc and E(B — V)—0A0 mag for the fore- 
ground population is shown by thick dashed curve. Fig. 5. 
reveals that the stars VI, V5, V14, V18, V21, V25, V29, 
V30, V35, V41 and V42 (shown with open circles) follow 
the ZAMS for field stars, hence could be MS foreground 
population. It is worthwhile to point that majority of the 
field stars (10 out of 11) are located outside the core (1.4 
arcmin, P08) of the cluster. 

In order to determine the ages and masses of the prob- 
able PMS variable stars we have used PMS isochrones and 
evolutionary tracks by Siess et al. (2000), respectively. Ma- 
jority of the probable PMS cluster members have masses in 
the range 0.3< M/Mq <3.5 and ages in the range l<age<5 
Myrs. This indicates that the young stellar objects asso- 
ciated with the cluster could be TTSs. Star V23 is the 
brightest variable in the sample and lies near 0.1 Myr PMS 
isochrone. 



3.2 NIR Colour-Colour Diagram 

The NIR colour-colour diagram (CCD) can be used to distin- 
guish the CTTSs and WTTSs. Fig. 6 shows NIR CCD for the 




Figure 5. V, V — I colour-magnitude diagram of Be 59. The 
ZAMS by Girardi et al. (2002) and PMS isochrones for 0.1, 1, 5, 
10 Myrs by Siess et al. (2000) arc also shown. The dashed curves 
show PMS evolutionary tracks for stars of different masses. The 
isochrones and evolutionary tracks are corrected for the cluster 
distance and E(B — V) = 1.45 mag. The filled circles and triangles 
represent probable WTTSs and CTTSs respectively (cf. section 
3.2). The thick dashed is ZAMS by Girardi et al. (2002) corrected 
for a distance of 470 pc and reddening E(B — V) = 0.40 mag. 
The open circles represent foreground stars. 
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Figure 6. (J — H, H — K) colour-colour diagram for stars lying 
in the 13x13 arcmin 2 field of Be 59. JHK data have been taken 
from 2MASS catalogue. The symbols are same as in Fig. 5. The 
sequences for dwarfs (solid) and giants (long dashed) are from 
Bessel & Brett (1988). The dotted line represents the intrinsic 
locus of CTTSs (Meyer et al. 1997). Dashed lines represent the 
reddening vectors (Cohen et al. 1981). 



identified variables. The NIR data have been taken from the 
2MASS Point Source Catalogues (PSC; Cutri et al. 2003). 
The 2MASS counterpart of the variables were searched us- 
ing a match value of 3 arcsec. Of the 42 detected vari- 
ables, we found 2MASS counterparts for 41 stars (30 cluster 
members and 11 foreground MS stars). The 2MASS coun- 
terparts for V32 star could not be found. It is the reddest 
variable star in the sample and lying above 0.1 Myr PMS 
isochrone in V, (V — I) CMD. However, on the basis of the 
location of this star on V, (V — I) CMD, we have considered 
V32 as a PMS object. The 2MASS magnitudes and colours 
were transformed to the CIT system using the relations 
given on the 2MASS website. The solid and long dashed 
lines represent unreddened MS and giant branch (Bessell & 
Brett 1988), respectively. The dotted line indicates intrin- 
sic locus of CTTSs (Meyer et al. 1997). The parallel dashed 
lines are the reddening vectors drawn from the tip (spectral 
type M4) of the giant branch ('upper reddening line'), from 
the base (spectral type AO) of the main-sequence branch 
('middle reddening line') and from the tip of the intrinsic 
CTTS line ('lower reddening line'). The extinction ratios 
Aj/A v = 0.265, Ah/Av = 0.155 and A K /A V = 0.090 have 
been adopted from Cohen et al. (1981). We classified the 
sources into three regions in the NIR CCD (cf. Ojha et al. 
2004a). The 'F' sources are those that are located between 
the upper and middle reddening lines and are considered to 
be either field stars (main-sequence stars, giants) or Class 
III and Class II sources with small NIR excesses. The iden- 
tification of Class II sources having a small amount of NIR 
excess is difficult on the basis of ( J — H, H — K) CCD. How- 
ever, recent studies, where Class II sources were identified on 
the basis of mid-IR observations through Spitzer telescope, 



manifest that these Class II objects lie in 'F' region near 
the middle reddening vector (e.g. Chauhan et al. 2011). The 
'T' sources are those that are located between the middle 
and lower reddening lines. These sources are considered to 
be mostly CTTSs (Class II objects). 'P' sources are those 
that are located in the region redward of the T region and 
are most likely Class I objects (protostar-like objects; Ojha 
et al. 2004b). 

Majority of the stars lying in the 'F' region are con- 
sidered as Class III stars/WTTSs and are shown by filled 
circles. However, there seems segregation of sources between 
(J-H) > 1.0 mag and (H — K) > 0.50 mag in {J—H, H-K) 
colour-colour diagram. As discussed above that Class II 
sources with small amount of NIR excess could be found 
in 'F' region near the middle reddening vector. The stars 
having (J — H) > 1.0 mag and (H — K) > 0.50 mag are con- 
sidered CTTS and these are shown by filled triangles. The 
foreground MS stars as identified from V, V — I CMD are 
distributed around the ZAMS and shown by open circles. 
The nature of the identified variables is mentioned in the 
Table 2. 

The membership of the variable stars is determined on 
the basis of their location in (V, V - 1) CMD and NIR CCD. 
However, spectroscopic follow up of these objects is needed 
to confirm their nature. 



4 PERIOD DETERMINATION 

We used the Lomb-Scargle (LS) periodogram (Lomb 1976; 
Scargle 1982) to determine the most likely period of a vari- 
able star. The LS method is useful for finding significant 
periodicities even with unevenly sampled data. We used 
the algorithm publicly available at the StarlinlQ software 
database, and verified the periods further with the software 
period04 (Lenz & Breger 2005). The software period04 
provides the frequency and semi-amplitude of the variabil- 
ity in a light curve. Periods derived from the LS method 
and Period04 generally matched well. In order to further 
verify the period estimates we also used periodogram analy- 
sis available onlinfl For any star showing a spurious period, 
we visually inspected the phased light curve for that particu- 
lar period. The phased light curves for CTTSs, WTTSs and 
field stars are shown in Figs. 7, 8 and 9, and 10 respectively. 
Figs. 7, 8, 9 and 10 plot the mean of differential magnitudes 
in 0.04 phase bin as a function of phase. The error bars rep- 
resent the a of the mean. In subsequent section we discuss 
the variability characteristics of individual variable star. 



5 VARIABILITY CHARACTERISTICS 
5.1 Non Periodic Variables 

Light curve of one WTTS variable named V38 displayed in 
Fig. 9 does not show periodic variability. This star could 
be either irregular variable or the period estimated in the 
present work could be in error. 

2 http://www.starlink.uk 

3 http:/ /www.univie.ac.at/tops/Period04 

4 http:/ /nsted. ipac.caltech.edu/periodogram/cgi- 
bin/Pcriodogram/nph-simpleupload 



Search for variable stars 



V3, p = 89.800, amp = 0.855 
-2 M 

O 

2 

O 0.5 1 1.5 2 

VIO, p = 58.753, amp=1.122 



O 0.5 1 1.5 2 

V17, p = 0.6B8, amp = 0.354 

o 
1 



J ' ' ' ' I 








%.f 




1 







O 0.5 1 1.5 2 

V34, p = 0.423, amp = 0.315 



-1 

o 
1 



J I 












1 





O 0.5 1 1.5 2 
Phase 



V4, p = 0.652 amp = 0.27'2 
-1 I | ' ' ' ' | ' ' ' ' | ■ ' ' ■ | ' ' ' ' 



Us 



o 
1 

O 0.5 1 1.5 2 

V13, p— 1.212, amp — 0.290 



2 


_U i i i | i i M | 
:■ . ■■ 




-0.5 


O 






O 


2 


T, , , , 1 , , , , 1 




0.5 









H , 







O 0.5 1 1.5 2 



V28, p = 9 735, amp = 0.356 
-1 I | ' ' ' ' | ' ' ' ' | ■ ' ' ■ | ' ' ' ' 



o 
1 



O 0.5 1 1.5 2 
Phase 



V7, p^6,586, amp = 0.3"70 



1st - " 



-0.5 

O 
0.5 

O 0.5 1 1.5 2 

V15, p=10,317 amp = 0.215 

-0.5 

O 
0.5 



T 1 



O 0.5 1 1.5 2 

V31, p = 0.516, amp = 0.221 

o 
1 



O 0.5 1 1.5 2 
Phase 



Figure 7. The phased light curves of probable CTTSs candidate variables identified in Be 59. 
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Figure 8. The phased light curves of probable WTTSs candidate variables identified in Be 59. 
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Figure 9. Same as Fig. 8. 
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Figure 10. The phased light curves of foreground MS variables. 
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5.2 Periodic Variables 

The short term periodic variations (2 ~ 10 days) in TTSs 
are believed to occur due to the axial rotation of a star 
with an inhomogeneous surface, having either hot or cool 
spots (Herbst et al. 1987, 1994). Herbst et al. (1994) iden- 
tified three types of TT variables based on their variability 
timescales ranging from a day to weeks. Type I variability, 
most often seen in WTTSs, is characterized by smaller stel- 
lar flux variations (a few times of 0.1 mag) and results from 
the rotation of a cool spotted photosphere. Type II varia- 
tions have larger brightness variations (up to ~2 mag), most 
often irregular but sometime periodic, are associated with 
short-lived accretion related hot spots at the stellar surface 
of CTTSs. The rare type III variations are characterized by 
luminosity dips lasting from a few days upto several months, 
which presumably result from circumstellar dust obscura- 
tion. The sample of TTSs in the cluster Be 59 indicates 
that ~90% TTSs have periods less than 15 days. The esti- 
mated periods of 8 CTTSs vary in the range of 0.652 day 
to 10.317 day, however stars V3 and V10 have significantly 
larger periods (89.8 day and 58.7 day respectively). The pe- 
riod estimates for WTTSs vary from 0.109 day to 5.671 day, 
however ~95% WTTSs have periods < 2 day. The shorter 
periods in the case of WTTSs could be explained as most of 
the stars at this age might have unlocked from their disks 
while contracting towards MS. They spin up because of de- 
creasing radii in order to conserve their angular momentum 
(Strom et al. 1989; Skrutskie et al. 1990; Haisch et al. 2001; 
Bouvier et al. 1994; Herbst & Mundt 2005; Lamm et al. 
2005; Kundurthy et al. 2006). But Nordhagen et al. (2006) 
have different view and suggest that environment also plays 
a crucial role in establishing the rotation periods besides 
age and mass. The presence of 'O' type stars in the vicinity 
would speed up the process of the dispersal of the circum- 
stellar disks causing the stars to spin up at a much younger 
age. 

The amplitude of WTTSs ranges from 0.022 to 0.273 
mag. The amplitude in the case of CTTSs has a range 
of 0.215 to 1.122 mag. The present results regarding am- 
plitudes are in agreement with Grankin et al. (2008 for 
WTTSs) and Grankin et al. (2007 for CTTSs). The bright- 
ness of CTTSs is found to vary with larger amplitude in 
comparison to the WTTSs. The large amplitude in the case 
of CTTSs indicates presence of the hot spots on their sur- 
faces. The smaller amplitude in the case of WTTSs suggests 
that these PMS stars might be in the process of dissipation 
of their circumstellar disk. 

The star V3, a probable CTTS, seems to be an inter- 
esting object. It has relatively high excess (H — K — 1.23) 
and shows relatively high extinction (cf. Fig. 6). It has a 
period of 89.80 days with amplitude variation of 1.0 mag 
and its mass is estimated to be ~0.9 Mq. The relatively 
high extinction of this object could be due to the presence 
of the circumstellar disk. Edwards et al. (1993) found that 
the stars which showed the evidence of circumstellar disks, 
inferred from their (H — K) colours, are slow rotators in 
comparison to those that lacked circumstellar disks. Rela- 
tively high extinction and longer period of V3 support this 
argument. This star may belong to type III as classified by 
Herbst et al. (1994), where relatively larger luminosity dips 



for longer period are expected due to circumstellar dust ob- 
scuration. 

The Ha emission star V10, classified as CTTS, also 
seems to be an interesting object. A careful inspection of 
night to night observations of this star reveals amplitude 
variation. Its 2006 observations show that the brightness of 
this star was varying with larger amplitude (1.122 mag) and 
period 1.102 day. From 2007 and onwards, it shows consis- 
tent brightness variation with amplitude of ~0.20 mag and 
period of 1.102 day. Whereas, the whole data set yield a 
period of 58.753 day. 

Star V19 ((V = 14.665 mag) could be a PMS pulsating 
star as suggested by its variability characteristics like period, 
amplitude and the shape of light curve. The period and am- 
plitude of star V19 vary over time. Its period on 11 October 
2009 and 27 October 2010 is found to be 0.101 day and 0.096 
day, respectively, whereas its amplitude varied from 0.022 to 
0.014 mag on respective dates. Considering all the observa- 
tions of this star its period comes out be 0.109 day. Star V19 
might be a 5 Scuti type variable star. 

An early type star V23 (brightest star in the present 
sample), considered to be a member of the cluster, might 
be a pulsating star. Its period (0.217 day) and amplitude 
(0.027 mag) are very much similar to the B-type pulsating 
stars. 

Star V21, lying just below the MS, could be a cluster 
member. The period, amplitude and the shape of light curve 
of V21 suggest that it might be a pulsating star. The spec- 
troscopic observations of star V21 was made by Majaess et 
al. (2008). They considered it to be a member of the cluster 
Be 59 and found it to a be B-type star. 

The field variables have period < 1.16 days and ampli- 
tude < 0.30 mag. 

5.3 Period and Amplitude Distribution 

In order to study the correlation between stellar activity and 
mass/age, we plot rotation period as a function of mass/age 
in Fig. 11. The rotation periods of TTSs associated with the 
cluster Be 59 do not show any dependency on the mass (Fig. 
11 left panel). Fig. 11 (right panel) shows two long periodic 
variables (V3 with P=89.8 day and V10 with P=58.753 day, 
cf. Table 2) have ages less than 5 Myrs, which could be inter- 
preted as the rotation is regulated by disk locking at early 
ages (Jayawardhana et al. 2006, Herbst et al. 2002). The 
presence of NIR excess in (H — K) of V3 and Ha emission 
in case of V10 is probably indicative of disk locking. But it 
is worth to mention here that large IR excess alone could 
not indicate ongoing disk-locking because IR excess is the 
indicator for the presence of circumstellar dust and it is not 
an indicator for current mass accretion (Lamm et al. 2005). 

Fig. 12 (left panel) and Fig. 12 (right panel) show am- 
plitude as function of mass and age, respectively. Fig. 12 
(left panel) reveals that the amplitude of TTSs seems to 
be correlated with the mass of the star in the sense that 
amplitude decreases with increasing mass of the variables. 
The largest amplitude is found in the case of V3 (0.855 mag) 
and V10 (1.122 mag). The star V3 shows significantly higher 
(H — K) excess with relatively larger extinction, whereas star 
V10 shows Ha emission. Thus both the stars manifest pres- 
ence of disk. The decrease in amplitude could be due to the 
dispersal of the disk material. If this notion is true, it in- 
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Figure 11. Rotation period as a function of mass and age of 
CTTSs and WTTSs. The symbols are same as in Fig. 5. The 
encircled triangles represent long and short periods for V10. 
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Figure 12. Amplitude as a function of mass and age of CTTSs 
and WTTSs. The symbols are same as in Fig. 5 and Fig. 11. Inset 
in the left panel shows a magnified view of amplitude variation. 



dicates that the mechanism for disk dispersal operates less 
efficiently for relatively low mass stars, fn fact there is ev- 
idence in literature that disk dispersal mechanism depends 
on the mass of the stars (e.g. Carpenter et al. 2006). But 
in Fig 12. (right panel), the amplitude variation with ages 
shows no definite trend, and majority of the sources having 
the ages within 5 Myrs show the amplitude of variability up 
to 0.35 mag with two CTTSs as exception, which have been 
discussed earlier. 



6 SUMMARY 

The paper presents our efforts to search for variable stars 
in young open cluster Be 59. Using VI photometry of a 
~ 13 x 13 field around the cluster we have identified 42 
variable stars. The probable members of the cluster are 
identified using (V, V — I) colour-magnitude diagram and 
(J — H , H — K) colour-colour diagram. Out of 42 variables, 
31 variable stars are found to be the probable PMS stars 
associated with the cluster Be 59. The majority of the prob- 
able PMS candidate variable stars have ages 1 to 5 Myrs. 
The masses of probable PMS variable stars range from ~0.3 
to ~3.5 Mq, suggesting that these could be TTSs. The sam- 
ple of TTSs in the cluster Be 59 indicates that ~90% TTSs 
have periods less than 15 days. The CTTSs have larger am- 
plitudes in comparison to WTTSs. The larger amplitude of 
CTTSs indicates the presence of hot spots on their surface. 
The smaller amplitude in the case of WTTSs suggests that 
these PMS stars might be in the process of dissipation of 
their circumstellar disk. It is found that the amplitude of 
TTSs decreases with increasing mass of the variable stars. 
The decrease in amplitude of variability of TTSs could be 
due to the dispersal of the disk material, indicating that 
the mechanism for disk dispersal operates less efficiently for 
relatively low mass stars. 
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